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Integrability of three dimensional models: cubic equations 
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We extend basic properties of two dimensional integrable models within the Algebraic Bethe 
Ansatz approach to 2+1 dimensions and formulate the sufficient conditions for the commutativity 
of transfer matrices of different spectral parameters, in analogy with Yang-Baxter or tetrahedron 
equations. The basic ingredient of our models is the R-matrix, which describes the scattering of a 
pair of particles over another pair of particles, the quark-anti-quark (meson) scattering on another 
quark-anti-quark state. We show that the Kitaev model belongs to this class of models and its 
R-matrix fulfills well-defined equations for integrability. 

PACS numbers: 


The importance of 2D integrable models^^® in mod¬ 
ern physics is hard to overestimate. Being initially an 
attractive tool in mathematical physics they became an 
important technique in low dimensional condensed mat¬ 
ter physics, capable to reveal non-perturbative aspects 
in many body systems with great potential of applica¬ 
tions. The basic constituent of 2D integrable systems is 
the commutativity of the evolution operators, the trans¬ 
fer matrices of the models of different spectral parame¬ 
ters. This property is equivalent to the existence of as 
many integrals of motion as number of degrees of freedom 
of the model. It appears, that commutativity of trans¬ 
fer matrices can be ensured by the Yang-Baxter (YB) 
equations'^”^ for the R-matrix and the integrability of the 
model is associated with the existence of the solution of 
YB-equations. 

Since the 80s of last century there was a natural desire 
to extend the idea of integrability to three dimensions®, 
which resulted in a formulation of the so-called tetra¬ 
hedron equation by Zamolodchikov^. The tetrahedron 
equations (ZTE) were studied and several solutions have 
been found until now^’®d 0 >i 3 ,i 4 ,i 7 ,i 9 ^ However, earlier so¬ 
lutions either contained negative Boltzmann weights or 
were slight deformations of models describing free parti¬ 
cles. Only in a recent work^® non-negative solutions of 
ZTE were obtained in a vertex formulation, and these 
matrices can be served as Boltzmann weights for a 3D 
solvable model with infinite number of discrete spins at¬ 
tached to the edges of the cubic lattice. In this sense it 
is remarkable to note that among the general solutions 
obtained in this paper it is also possible to detect R- 
matrices with real and non-negative entries which can be 
considered as Boltzmann weights in the context of the 
3D statistical solvable models with 1/2-spins attached to 
the vertexes of 3D cubic lattice. 

Although initially the tetrahedron equations were for¬ 
mulated for the scattering matrix S of three infinitely 
long straight strings in a context of 3D integrability they 
can also be regarded as weight functions for statistical 


models. In a Bethe Ansatz formulation of 3D mod¬ 
els their 2D transfer matrices of the quantum states on 
a plane®’can be constructed via three particle R- 
matrix®’^'*’^' , which, as an operator, acts on a tensorial 
cube of linear space V, i.e. R-.V®V®V ^ . 

Another approach to 3D integrability based on Frenkel- 
Moore simplex equations^® also uses three-state R- 
matrices. They are higher dimensional extension of quan¬ 
tum Yang-Baxter equations without spectral parameters. 
However these equations are less examined^"^. 

Motivated by the desire to extend the integrability con¬ 
ditions in 3D to other formulations we consider a new 
kind of equations with the R-matrices acting on a quar- 
tic tensorial power of linear spaces V 

.Ri 234 :Fl(8>V2®V3®V4—>'V"i(8>V2®V3® V4, (0.1) 

which can be represented graphically as in Fig. la. The 
R-matrix can be represented also in the form displayed 
in Fig. lb, where the final spaces are permuted (Vi and 
V 2 with V 3 and V4, respectively): R1234 = ^I234^’i3^24- 
Explicitly it can be written as follows 

d/3i/32/33/34 _ d/33/34/3i/32 /q o') 

■*^^0:10:20:30:4 ^0:10:2030:4 * V ^’^) 

Identifying the space Vi ( 8 > V 2 and V 3 ( 8 14 with the quan¬ 
tum spaces of quark-anti-quark pairs connected with a 
string one can regard this R-matrix as a transfer matrix 
for a pair of scattering mesons. Within a terminology 
used in the algebraic Bethe Ansatz for 1+1 integrable 
models this R-matrix can be viewed also as a matrix, 
which has two quantum states and two auxiliary states. 

The space of quantum states of 

the system on a plane is defined by a direct product of 
linear spaces of quantum states on each site (n, m) 
of the lattice C (see Fig. 2a). We fix periodic bound¬ 
ary conditions on both directions: 14 ,,m+L = Vn,m and 
Vn+L,m = Vn,m- The time evolution of this state is deter¬ 
mined by the action of the operator/transfer matrix T: 
= ^tT, which is a product of local evolution opera- 
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FIG. 1: Four particle R-matrix 


tors, R-matrices as follows. First we fix a chess like struc¬ 
ture of squares on a lattice C and associate to each of the 
black squares a R-matrix R(„+i,„)(„+!,„+!)(„,„)(„,„+i), 
which acts on a product of four spaces at the sites. In 
this way the whole transfer matrix becomes 


T = TrU: 


LI 2 


n—1 
L/2 n 


jjL/2 p 

-*--^m=l'^(2n,2m){2n,2m+l)(2n—l,2m)(2n—l,2m+l) 


• ^rn=l^i2n+l,2m-l)(2n+l,2m){2n,2m-l){2n,2m) ; (0.3) 


where the Trace is taken over states on boundaries. 
The indices of the R-matrices in the first and sec¬ 
ond lines of this product just ensure chess like order¬ 
ing of their action. In Fig. 2b we present this prod¬ 
uct graphically. First we identify the second pair of 
states {{2n — 1, 2m), { 2 n — 1, 2 m -I- I)) (in first row) and 
((2n-|-I, 2m—l)(2n-|-I, 2m))(in second row) of R-matrices 
with the corresponding links on the lattice. Then we ro¬ 
tate the box of the R-matrix by 7r/4 in order to ensure the 
correct order for their action in a product. In the same 
way we define the second list of the transfer matrix, which 
will act in the order TbTa- Fig. 2c presents a vertical 2D 
cut of two lists of the product TbTa drawn from the side. 
The 7r/4 rotated lines mark the spaces 14,m attached to 
sites (n, m) of the lattice. Though transfer matrix (0.3) 
is written in R formalism, it can easily be converted to 
the product of i?-matrices. 

The arrangement of R-matrices in the first row (first 
plane of the transfer matrix Tg) acts on the sites of dark 
squares of the lattice while R-matrices in the second row 
(second plane of the transfer matrix Ta) act on the sites 
of the white squares. 

Being an evolution operator the transfer matrix 
should be linked to time. According to the general 
prescription"*’^ the transfer matrix T{u) is a function 
of the so-called spectral parameter u and the linear 
term Hi in its expansion T(u) = defines the 

Hamiltonian of the model, while the partition function 
is Z = TrT^. Integrable models should have as many 
integrals of motion, as degrees of freedom. This property 
may be reached by considering two planes of transfer ma¬ 
trices with different spectral parameters, T(u) and T(v) 
and demanding their commutativity [T{u)^T{v)] = 0, or 
equivalently demanding the commutativity of the coeffi¬ 
cients [Hr, Hs] = 0 of the expansion. This means, that all 





(n.m) 


FIG. 2: Transfer matrix T as a product of R-matrices on 
the plane: a) presents the lattice with chess like structure; 
b) presents the product of R-matrices arranged in chess like 
way; c) 2D cut of the product TbTa- 


Hr, r > 1 are integrals of motion. In 2D integrable mod¬ 
els the sufficient conditions for commutativity of trans¬ 
fer matrices are determined by the corresponding YB- 
equations'^“®. 

In order to obtain the analog of the YB equations, 
which will ensure the commutativity of transfer matrices 
(0.3) we use the so-called railway construction. Let us 
cut horizontally two planes of the R-matrix product of 
two transfer matrices (on Fig. 2b we present a product of 
R-matrices for one transfer matrix plane) into two parts 
and substitute in between the identity 


n 


L 

m—1 


TrU 


Tr 


n 


^^(2n-t-l,m)^^(2n,Tn) 

m—l"^(2n-t-l,m)(2n-t-l,m-t-l)(2n,m)(2n,m-t-l) 
m—1 -^(2n-t-l,m) (271-1-1,771-t-1) (2n,m)(2n,m-t-l) 


(0.4) 


which maps two chains of sites, {2n,m),m = \---L 
and (2n -|- l,m),m = 1 ■ ■ ■ L + \, into itself. The 
Trace have to be taken by identifying spaces I and 
L-l-I. In this expression we have introduced another 
set of i?-matrices, called intertwiners, which will be 
specified below. For further convenience we distinguish 
-^( 2 n-t-i,m)( 2 n-t-i,m-t-i)( 2 n,m)( 27 i, 77 i-t-i) matrices for even and 
odd values of m marking them as i ?3 and i ?4 respectively. 
In the left side of Fig. 3 we present one half of the plane 
of R-matrices together with an inserted chain of R 3 R 4 
as intertwiners. The chain of intertwiners can also be 
written by i?-matrices. 

Now let us suggest, that the product of these inter¬ 
twiners with the first double chain of R-matrices from 
the product of two planes of transfer matrices is equal 
to the product of the same operators written in opposite 
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FIG. 3: Reduced set of commutativity conditions for transfer 
matrices. The dotted subset represents the cubic equations 
( 0 . 6 ). 

order. Namely we demand, that 

him—l'^(2n+l,m) (2n+l,m+l) (271,771) (2n,m+l) 

• n,^'^j^i?(277,2m)(277,2m+l)(277-l,2m)(27i-l,2m+l)(u) 

■ h[^'^j^.R(277+1.2m)(277+l,2m+l)(277,2m)(27i,2m+l) (v) 

= ^m=l^{2n,2m){2n,2m+l){2n-l,2m){2n-l,2m+l)i'>^) 

■ h[,^'^j^i?(277+1.2m)(27i+l,2m+l)(277,2m)(277,2m+l) (w) 

■^771^-^(271-1-1 ,?n) (2?T-1- 1,771-1-1 ) ( 271 , 777 ) (271,771-1-1) ’ (0.5) 

Graphically this equation is depicted in Fig. 3. We move 
the column of intertwiners from the left to the right hand 
side of the column of two slices of the R-matrix product, 
simultaneously changing their order in a column, chang¬ 
ing the order of spectral parameters u and v of the slices 
and demanding their equality. We can use the same type 
of equality and move the chain of intertwiners further to 
the right hand side of the next column of the two slices 
of the i?-matrix product. Then, repeating this operation 
multiple times, one will approach the chain of inserted 
R~^ intertwiners inside the Trace from the other side 
and cancel it. As a result we obtain the product of two 
transfer matrices in a reversed order of spectral param¬ 
eters u and V. Hence, the set of equations (0.5) ensures 


the commutativity of transfer matrices. 

The set of equations (0.5) can be simplified. Namely, 
it is easy to see, that the equality can be reduced to the 
product of only 2 7?-matrices, R{u) and R{v) and two 
intertwiners, R 3 and R 4 . In other words, it is enough 
to write the equality of the product of 7?-matrices from 
the inside of the dotted line in Fig. 3. Graphically this 
equation is depicted in Fig. 4. 

We see, that in this equation the product of R-matrices 
acting on a space (for simplicity we numerate the 

spaces from 1 to 9) can be written as 

'^5263('*^J ^)'^4152('^J ^)^5689('^)'^2356(^)*'^7 

= -^4578('^)-^1245('*^)^l596('*^i'^)-^7485('*^i(0-6) 

Here we have introduced a short-hand notation for R- 
matrices simply by marking the numbers of linear spaces 
of states, in which they are acting; id^ and id^ are identity 
operators acting on spaces 3 and 7 respectively. Eq. (0.6) 
can also easily be written by use of R. 

This is the set of equations, sufficient for commutativ¬ 
ity of transfer matrices. The same set of equations are 
sufficient for commuting i?-matrices in the second column 
in Fig.l. Equations (0.6) form an analog of YB equations 
ensuring the integrability of 3D quantum models. Since 
they have a form of relations between the cubes of the 
R-matrix picture(see Fig.l) we call them cubic equations. 

We will show now that the Kitaev inodel^^ can be de¬ 
scribed as a model of the prescribed type and its R- 
matrix fulfills the set of cubic equations (0.6). The 
full transfer matrix of Kitaev model is a product TaTb 
of two transfer matrices of type (0.3) defined by R- 
matrices Ra = Ox ® <Jx ® <yx ® 

and Rb = \ + u (Jz®(yz®crz®(Jz respectively. 

The linear term of the expansion of TaTb in the spectral 
parameter u will produce the Kitaev model Hamiltonian 

^Kitaev — ^ ( ^x ® ^x ® 

white plaquettes 

+ ^ CTz ® CTz ® (Jz ® <Jz- (0.7) 

dark plaquettes 


The integrability of the Kitaev model is trivially clear 
from the very beginning since all terms in the Hamilto¬ 
nian defined on white and dark plaquettes commute with 
each other. The latter indicates, that the number of in¬ 
tegrals of motion of the model coincides with its degrees 
of freedom. However, in this paper we are aiming to 
show, that one can develop 3D Algebraic Bethe Ansatz 
approach in such a way, that the Kitaev model will auto¬ 
matically be integrable. Namely, we will show now, that 
Ra and Rs-matrices of the Kitaev’s model fulfill Eq. 0.6. 
The explicit form of Eq. 0.6 by use of indices according 


to the definition in Fig. la reads 

537471^5/34 52^’’^*7879 , N 5372737576 , , e77 _ 

^ a^azagaa^ a4ai/3i/32'^ ^ PiPgPzPg\^>^az 

-^PiPsPaPa . .-^in^PiPs , , -3^1^27976 n47'’''7^'7®'75 73,3 
^ 040:50708 0 1 02/37/38 1^1"^ PgPgagag^ P^PiPlP^^O-g- 

(0.8) 

where R^{u) = Ra{u) and R'^{v) = Rb{v)- It appears, 
that the intertwiners 

i?4 = Ra~\u), = RBiv) (0.9) 

where R']^{u) = I®!®!®! — uux ® cfx ® Ux ® cfx 
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FIG. 4: The set of equations ensuring commntativity of transfer matrices with different spectral parameters. We numerate 
linear spaces of states V, where R-matrices are acting, by Vi • • • Vg. R-matrices on the left hand side of equation are not acting 
on space V 7 while R-matrices on the right hand side are not acting on V 3 . We pnt identity operators in the equation acting on 
this spaces for the consistency. 


fulfill the cubic equations (0.8) for any parameters u and 
V. This can be directly checked both, by a computer 
algebra program and analytically. The commutativity 
of transfer matrices Ta(u) with Ta{v) and Tb{u) with 
Tb{v) is trivial in the Kitaev model. 

Summary. We have formulated a class of three di¬ 
mensional models defined by the i?-matrix of the scat¬ 
tering of a two particle state on another two particle 
state, i.e. a meson-meson type scattering. We derived 
a set of equations for these i?-matrices, which are a suf¬ 
ficient conditions for the commutativity of the transfer 
matrices with different spectral parameters. These equa¬ 


tions differ from the tetrahedron equations, which also 
ensure the integrability of 3D models, but are based on 
the R-matrix of 3 particle scatterings. Our set of equa¬ 
tions will be reduced to tetrahedron type of equations 
by considering the two auxiliary spaces in the R-matrix 
as one (fusion) and replacing it by one thick line. We 
showed that the Kitaev modeD^ belongs to this class of 
integrable models. 
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